Mon. Not. R. Astron. Soc. 000, 1-?? () Printed 5 February 2008 (MN J^I^X style file v2.2) 



Cold Dust in (Some) High z Supernova Host Galaxies 

D.L. Clements 1 , D. Farrah 2 , M. Rowan-Robinson 1 , J. Afonso 3 ,R. Priddey 4 , M. Fox 1 

1 Imperial College, London, Blackett Lab, Prince Consort Road, London SW7 2BW, UK 
2 IPAC, Caltech, 770 S. Wilson Ave., Mailstop 100-22, Pasadena, CA 91125, USA 

3 CAAUL, Universidade de Lisboa, Faculdade de Ciencias,Observat6rio Astronomico de Lisboa, Tapada da Ajuda, 1349-018 Lisboa 
4 Department of Physics, Astronomy and Mathematics, University of Hertfordshire, Hatfield, Herts., UK 



ABSTRACT 

We present deep submillimetre photometry for 14 galaxies at z=0.5 that are hosts 
of type la supernovae, with the aim of examining the evolution of dust mass and 
extinction in normal galaxies. We combine these results with our previous observations 
of 17 z^0.5 SNla hosts to look for any evolution in the dust content of normal galaxies 
between z=0 and z=0.5. The average observed frame 850/im flux of SNla hosts in the 
full sample, excluding 2 bright individually detected objects, is 0.44 ± 0.22 mJy. This 
flux level is consistent with there being little or no evolution in the dust content, 
or optical extinction, of normal galaxies from z=0 to z=0.5. One galaxy, the host of 
SN1996cf, is detected individually, and we also present a deep HST STIS image for 
this object. It appears to be an edge on disk system, similar to the submm bright host 
of SN1997ey. We thus examine the dust properties of these and one other individually 
detected object. 450-to-850 ^m flux ratios and limits suggest that the dust in the 
two brightest submm sources, SN1996cf and SN1997ey, is cold, T ^20K, implying 
that they contain a substantial mass of dust ~ 10 9 M Q . The presence of two bright 
(Fsso >7mJy) submm sources at z~0.5 in a sample of ostensibly normal galaxies is 
surprising, and has important implications. It supports the idea that a substantial 
part of the Cosmic Infrared Background (CIB) may be produced at z< 1, while also 
suggesting that 'foreground' objects such as these may be a significant 'contaminant' 
in submm surveys. Finally, we examine the overall submm luminosity distribution at 
z=0.5 implied by our results, and conclude that either there is substantial evolution 
in the submm luminosity function from z=0 to 0.5, or our submm detected sources 
are somehow not representative of the bulk of galaxies at this redshift. 
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1 INTRODUCTION 

Two of the major new insights of the 1990s into cosmology 
and the history of galaxy formation were the discovery of 
the Cosmic Infrared Background (CIB: Puget et al., 1996; 
Fixsen et al., 1998) and the measurement of a significant 
'dark energy' term in the expansion of the universe through 
the use of high redshift type la supernovae as standard can- 
dles (Riess et al., 1998, HZT; Perlmutter et al., 1999, SCP). 
The discovery of the CIB indicated that dust enshrouded 
star formation is a significant aspect of the star formation 
history of the universe, amounting to 50% or more of all star 
formation (Gispert et al., 2000). At the same time surveys 
with SCUBA (eg. Smail et al., 1997; Hughes et al., 1998; 
Eales et al., 2000; Mortier et al., 2005) and other submillime- 
tre array detectors have begun to find the objects that make 
up the CIB. These objects are largely interpreted as being 
similar to local Ultraluminous Infrared Galaxies (ULIRGs), 
and are thought to contain large dust masses at a temper- 



ature T~40K (Blain et al., 1999) and to be the hosts of 
massive bursts of star formation, 100 — 1000Mot/r -1 . How- 
ever, substantial uncertainties remain. Most of these sub- 
millimetre galaxies (SMGs) do not have well determined 
redshifts (though see Chapman et al., 2005). The degen- 
eracy between temperature and redshift (Blain 1999) thus 
means that their dust temperature is highly uncertain — 
a low dust temperature object at low redshift (z~0.5 — 1) 
looks much the same as a high dust temperature object at 
high redshift (z~ 2 — 3). Since luminosity is a strong power 
of temperature, typically T 6 for a standard SMG spectral 
energy distribution (SED), this leads to a considerable un- 
certainty in the derived luminosity and star formation rates 
for those SMGs without measured redshifts. Indeed some 
authors have suggested (Rowan-Robinson 2001, Kaviani et 
al. 2003, Efstathiou & Rowan- Robinson, 2003) that a sub- 
stantial fraction of the SMG population may be cooler and 
closer than originally thought, and there is some observa- 
tional evidence to back up this suggestion (Chapman et al., 
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2003; Taylor et al in prep). There is currently little known 
about cool dust in normal galaxies in the 0.5<z<l range 
because it is quite difneuly to find a clean sample of such ob- 
jects. The SNla host galaxies (see eg. Sullivan et al. (2003); 
Tonry et al. (2003)), though, provide an ideal selection for 
such studies. 

The discovery of the CIB also raises the overall ques- 
tion of the evolution of the dust content of galaxies and 
its role in obscuring starlight in these objects and repro- 
cessing it into thermal dust emission. Dust obscuration is 
known to have a significant effect locally (Tresse & Maddox, 
1998), but the evolution of dust content and obscuration is 
not well constrained. Chemical evolution models (Calzetti 
& Heckman, 1999; Pei, Fall & Hauser, 1999) predict that 
dust obscuration should peak at 2-3 times the current value 
at 1 < z < 2 before gradually declining to higher redshifts. 
However, we currently know little about the dust content 
of quiescent systems which make up the bulk of the galaxy 
population at any epoch. The host galaxies of high z su- 
pernovae, however, have not been selected on the basis of 
active star formation, but have spectroscopic redshifts and 
abundant ancillary data, including, for many, HST images 
and multicolour optical and near infrared colours. They are 
thus an ideal sample with which to study the role of dust in 
quiescent galaxies at cosmological redshifts. 

This paper is the second to result from our ongoing 
search for dust emission in high z SN host galaxies. Our 
first paper (Farrah et al., 2004, hereinafter Paper 1) pre- 
sented the results of SCUBA submm photometry of a sample 
of 17 z=0.5 SNla host galaxies. We here present additional 
observations of a further 14 galaxies, extending the sample 
size to 31, and increasing the number of directly detected 
objects to three. The paper is structured as follows. In sec- 
tion 2 we discuss the observations and present our results. 
In section 3 we discuss the sample as a whole and discuss 
the overall dust content of the quiescent galaxy population 
at z=0.5. In section 4 we examine the properties of the three 
galaxies individually detected in our observations, while in 
section 5 we discuss the implications of our results for the 
overall submillimetre luminosity distribution of galaxies at 
z=0.5. We draw our conclusions in section 6. We assume 
fio = l,fi m = 0.3, «a = 0.7 and H = IQkms^Mpc- 1 
throughout this paper. 



2 OBSERVATIONS 

The observations described here extend the earlier work of 
Paper 1 on the host galaxies of z~0.5 type la supernova host 
galaxies. As with our Paper 1 observations, our targets were 
selected from the HZT and the SCP supernova lists to be the 
host galaxies of type la supernovae and to lie in a narrow 
(± 0.1) range of redshifts centred on z=0.5. This is to allow 
easy combination of submm fluxes from the targets with no 
need for fc-corrections to shift the fluxes to correspond to 
the same rest-frame wavelength. 

Observations were made of 14 high redshift supernova 
la host galaxies at z~ 0.5 between July 2003 and January 
2004 (see table 1). Observing conditions where generally 
good, with T850 ranging from 0.2 to 0.15 and T450 ranging 
from 1.3 to 0.85. These optical depth values were determined 
by standard skydip observations at azimuths relevant to the 



targets and are classified as grade 1 or 2 conditions by the 
JCMT. 

The sources were observed with SCUBA (Holland et 
al., 1999) in photometer mode using two bolometer chop- 
ping. This means that for each channel, 850 or 450 jj,m, 
as well as being observed by the central bolometer on the 
array, a second bolometer observes the source in the refer- 
ence position. An observed flux can be extracted from both 
these bolometers which are then combined, using appropri- 
ate noise weighting, to increase the sensitivity over that of 
a normal single bolometer chopping observation by ~ y/2. 
This combination takes place at the end of the data reduc- 
tion process, once flat fielding, background subtraction and 
calibration have been completed. 

The data from SCUBA was reduced in a standard way 
using the SURF data reduction package (Jenness & Light- 
foot, 2000). Regular pointing observations were made to en- 
sure pointing accuracy, and skydips were also taken regu- 
larly. Flux calibration factors (FCFs) for the two bolome- 
ters used to detect the sources were calculated separately. 
Calibration sources used were CRL618 and Uranus. The in- 
dividual observations of a source with a given bolometer 
are combined using a Kolmogorov-Smirnoff technique which 
masks out any discordant points that remain after the data 
processing. No unusual numbers of points were excluded in 
this analysis. When an object was observed on more than 
one day, the fluxes measured in each separate observation 
are combined using a variance weighted scheme. 

Our results are presented in table 1. Only one source, 
the host galaxy of SN1996cf, a z=0.57 SNla, was detected, 
with a flux of 11±1.6 mjy at 850pim. This source was only 
'detected' with 1.6<r significance at 450/im, with a flux of 
25±15 mjy. 

For the host galaxy of SN1996cf, we obtained Hubble 
Space Telescope (HST) imaging from the HST data archive. 
Observations were taken using the Space Telescope Imaging 
Spectrograph (STIS) in imaging mode using a clear filter. 
The data consisted of three exposures, each of 300 seconds, 
taken with a 10 pixel offset between successive exposures to 
facilitate the subtraction of cosmic ray events. The data were 
reduced, and combined into a single image, using the IRAF 
reduction package calstis. The image of the host galaxy, pre- 
sented in Fig. 1, has had the low surface brightness host 
enhanced by the application of a 2x2 pixel boxcar filter. 



3 THE EVOLUTION OF OBSCURATION IN 
NORMAL GALAXIES 

The evolution of dust content and dust obscuration in nor- 
mal galaxies can be examined by looking at the average 
properties of the sample we have observed. We follow the 
approach discussed in Paper 1, but apply it to the full sam- 
ple of 31 galaxies — 17 from Paper 1 and 14 from the cur- 
rent paper. We exclude 2 objects from this combination, the 
two objects whose fluxes are strongly detected — SN1996cf 
(this paper) and SN1997ey (Paper 1). The fluxes are com- 
bined using an inverse variance weighting scheme. The aver- 
age observed frame 850pim flux from type la supernova host 
galaxies calculated in this way is 0.44 ± 0.22 mjy. This is 
less than the mean flux of 1.01±0.33 mjy obtained in Paper 
1, but still statistically consistent (at the 1.5a level) with the 
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Figure 1. STIS clear CCD image of the host galaxy of SN1996cf. 
. The pixel size in this the standard STIS CCD resolution of 0.05". The Image is 15"xll" in size, ic. 300x220 pixels. 



Name 


z 


Date Observed 


Fs5o(mJy) 


SN1995az 


0.45 


9 July 2003 


1.4±0.9 


SN1997ax 


0.615 


26 Dec 2003 


1.8±0.8 


SN1997K 


0.592 


16 Dec 2003, 9 Jan 2004 


-0.3 ±0.9 


SN1997H 


0.526 


8, 10 Jan 2004 


0.7 ±0.9 


SN1997eq 


0.538 


8 Jan 2004 


1.5±1.2 


SN1997j 


0.619 


8 Jan 2004 


0.4 ±0.9 


SN1997af 


0.579 


9, 10 Jan 2004 


0.0±0.8 


SN2000cc 


0.47 


10 Jan 2004 


-2.7±1.4 


SN1996au 


0.52 


10 Jan 2004 


-3.0±1.4 


SN1997es 


0.65 


10 Jan 2004 


-0.7±1.4 


SN1996cf* 


0.57 


10 Jan 2004 


11.0±1.6 


SN1996ci 


0.5 


11, 12 Jan 2004 


-0.7±0.8 


SN1997aj 


0.581 


12 Jan 2004 


-1.8±1.3 


SN1996I 


0.57 


12 Jan 2004 


-0.6±1.2 



Table 1. Results of submm observations of SNla host galaxies. 
*SN1996cf 450Atm flux 25±15 mjy. 



previous result. If we attempt to split the SN host galaxy 
sample into spirals (and irregulars) and ellipticals, based on 
the available data on the host galaxy morphologies (largely 
HST images), whilst still excluding the two strong detec- 
tions, we find a similar average flux for the two classes of 
object (0.40 ±0.29 for spirals/irregulars, and 0.69 ± 0.46 for 
ellipticals) . 

By consideration of the local SLUGs 850/im luminosity 
function (Dunne et al., 2000 - though see also Vlahakis et 
al., 2005) combined with appropriate k corrections from the 
cirrus model of Efstathiou & Rowan- Robinson (2003), and 
the assumption that the submillimetre emission from ellipti- 
cals is negligible, Paper 1 predicted that the mean observed 



frame 850/xm flux of a z=0.5 galaxy would be 0.56 ±0.1 
mjy if there is no evolution in the amount of dust in normal 
galaxies from z=0 to z=0.5. Our result is thus consistent 
with there being at most only moderate evolution in the 
dust content of normal galaxies from z=0 to z=0.5. Paper 1 
found a moderate increase in the dust content of 25%-135% 
over this redshift range. The new expanded SN host sample 
thus mildly contradicts this result (1.5 a), though the errors 
permit an increase of up to 20% or a 100% decrease in the 
dust content of the population as a whole. The mean Ay for 
galaxies at z=0.5 would thus be statistically unchanged from 
the local value of Av = 0.31 derived by Rowan-Robinson 
(2003). This is in contrast with models of opacity evolution 
from Calzetti & Heckman (1999), which predict an increase 
in extinction of Av ~ 0.15 by z=0.15, and with chemical 
evolution models (Pei et al., 1999). Neither of these models 
predict a particularly significant evolution in extinction or 
dust mass in galaxies over the period to z=0.5, so our cur- 
rent results cannot be interpreted as a rigorous test of these 
models. To achieve this, a larger sample and more sensitive 
observations would be needed, allowing more than statis- 
tical non-detections of the population, and/or the greater 
lever arm that would be obtained by examining the dust 
content of normal galaxes out to z=l. 



All of these considerations apply to the population as a 
whole. Throughout this discussion we have excluded the two 
strong submm detections in our sample from consideration 
because the two objects appear to have radically different 
submm luminosity to the bulk of the population. We now 
consider the nature of these two objects in detail. 
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4 THE PROPERTIES OF THE DETECTED 
GALAXIES 

Whilst we have failed to detect the mean submm emission of 
the population of z~0.5 SNla host galaxies, we have man- 
aged to detect three of these objects individually. The SN 
magnitudes for these objects lie in typical positions for z~0.5 
SNla objects in the magnitude-redshift diagram (eg. Perl- 
mutter et al., 1999). One of our sources, SN2000eh, was 
only marginally detected (3u) at 850/im, but two others, 
SN1997ey and SN1996cf were detected surprisingly strongly. 
Their observational and derived properties are given in Table 
2. It can clearly be seen that their rest-frame 850/im lumi- 
nosities are comparable with those of local ULIRGs such as 
Arp220 and Mrk231 (Dunne et al., 2000; Farrah et al., 2003). 
However, there is no evidence for starburst or AGN activity 
in these objects on the basis of optical spectra obtained by 
the high redshift supernova teams (HZT & SCP teams, pri- 
vate communication). HST imaging of the two bright submm 
sources also shows them to be morphologically different from 
local ULIRGs and from those SMGs for which imaging is 
currently available. Both of our bright submm sources (see 
fig. 1 for 1996cf and Paper 1 for 1997ey) appear to be faint 
quiescent edge-on disk galaxies. Local ULIRGs, in contrast, 
are almost universally disturbed systems (eg. Clements et 
al., 1996; Surace et al., 1998; Farrah et al., 2001), with dou- 
ble nuclei, tidal tails and other signs of merging activity. 
The SMG population, as revealed by blank field SCUBA 
surveys, also appears to be made up of such disturbed ob- 
jects (eg. Clements et al 2004; Smail et al., 2004; Chapman 
et al., 2003b), though it should be noted that to date no flux 
limited SMG sample has been completely imaged. 

The spectral energy distributions and dust tempera- 
tures of these objects are difficult to determine since they 
are detected in only one or two submm bands. The best we 
can attempt is to set constraints on the temperature and 
SEDs by examining the 450/850 flux ratio and its limits. 
These targets were all observed in excellent conditions, so 
calibration errors should be smaller than or comparable to 
the observational errors on the fluxes (see eg. Dunne & Eales, 
2001). One source is clearly detected at 450/xm, SN1997ey, 
while a second, 1996cf, has a marginal 1.6 a signal. There 
is no useful 450/im data available for SN2000eh since it was 
observed in poorer conditions. Examination of the 450/850 
flux ratio for 1997ey and 1996cf (see fig. 2) suggest cold 
dust temperatures ~20K assuming an emissivity index of 
/3=1.3 (consistent with the single temperature SEDs found 
for local galaxies by the SLUGS survey (Dunne et al, 2000)) 
and still colder for the more generally accepted emissivity 
index (3=2 (see eg. Dunne & Eales 2001). Assuming a stan- 
dard dust opacity of 0.077 m 2 kg _1 (eg. Dunne et al. 2000), 
these would correspond to dust masses of 1.3xlO 9 M and 
1.7 xl0 9 M Q respectively. In the absence of 450/im data, we 
assume that SN2000eh has similar dust properties, yielding 
a lower dust mass of 6.5 x 10 Mq . These are substantial dust 
masses, about a factor of 10 more dust than in the dusti- 
est galaxy in the original SLUGS survey, UGC9618 (Dunne 
et al., 2000). The low inferred temperature, combined with 
the large observed 850/im flux, lead inevitably to this large 
dust mass. If the cooler temperature associated with a (3 = 2 
SED is assumed, then the dust masses get correspondingly 
larger. 



Dust as cold as we have found in these objects is not 
unknown in the local universe (eg. the SED of NGC891 
is dominated by dust at ~ 15K (Alton et al., 1998)) and 
has been suggested for higher redshift objects on both 
theoretical grounds (Rowan-Robinson, 2001; Efstathiou & 
Rowan-Robinson 2003; Kaviani et al., 2003) and observa- 
tional grounds (Chapman et al., 2003; Taylor et al., 2005). 
The present results are direct confirmation that cold galaxies 
with large dust masses and submm luminosities do indeed 
exist at moderate redshift, and could contribute to the CIB 
and to the sources seen in deep far-IR surveys such FIR- 
BACK (Puget et al., 1999) and deep submm surveys. Indeed, 
objects such as these might prove to be a foreground contam- 
inant in deep large area submm surveys such as SHADES 
(Mortier et al., 2005) and those now been planned for new 
facilties and instruments such as SCUBA2 and the LMT 
(Large Millimetre Telescope). The areal density for sources 
such as the two bright galaxies discussed here was predicted 
in Paper 1 to be ~100 per sq. deg. The results of this paper, 
which roughly doubles the sample size and doubles the num- 
ber of bright sources found, adds credence to this number, 
which is comparable to the source densities found in deep 
blank field submm surveys at these flux levels (eg. Scott et 
al., 2002; Mortier et al., 2005). The presence of such fore- 
ground interlopers might thus be a significant problem when 
attempting to measure correlation functions for high redshift 
submm galaxies in the absence of clear identifications and 
redshifts, especially as they are unlikely to be bright in the 
radio. 



5 Z=0.5 SUBMM LUMINOSITY 
DISTRIBUTION 

In principle, our sample of SNla host galaxies should pro- 
vide a reasonably unbiased selection for evolved galaxies 
at z=0.5, biased in selection only by the stellar mass. By 
examining the 850/im luminosity distribution of these ob- 
jects, and comparing it with that found for local objects, we 
can gain some insight into the galaxy population at z=0.5. 
The local luminosity function that we use is the SLUGS 
850/mi LF (Dunne et al 2000) since this is the best currently 
available. This was found to be well fitted by a Schechter 
function with L„ = 8.3 x 10 21 IUm~ 2 sr _1 , a = -2.18 and 
$» = 2.9 x 10~ 4 M pc~~ 3 dex^ 1 We cannot apply this luminos- 
ity function directly to predict the numbers of sources we 
should see, since we do not know the parent population nor- 
malisation for the number of our target objects per Mpc~ 3 . 
However, since all our objects are at the same redshift, ~0.5, 
we can look at the relative number of objects in each lumi- 
nosity bin and compare this to the low z expectation. We 
have only a small number of detected galaxies, so must nec- 
essarily normalise the expected luminosity distribution at 
the high end, where our detections lie. 

The Schechter function form of the LF used by Dunne 
et al. (2000) has a very strong exponential rolloff at high 
luminosities. The detection of a small number of very lumi- 
nous objects, as obtained here, would imply the presence of 
many more lower luminosity systems. We find 3 sources in 
the luminosity range 10 22 ' 5 - 10 23 ' 5 Wm^sr^ 1 . The Dunne 
et al. LF would predict that the next lowest dex bin in the 
LF should contain 100 times more objects. These would have 
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Figure 2. Comparison of Dust SED models to 450/850 ratio of strongly detected sources 

The curved lines are the predicted observed-frame 450/850 ratios for thermal dust emission models with a /3=1.3 emissivity (dashed 
line) or j3=2 emissivity (solid line). The horizontal lines are the observed ratios. Solid lines show the observed values and, for SN1996cf, 

the upper limit to this ratio using the formal 450/im 3<r flux upper limit (designated with an arrow), while dashed lies show the la 
range for this flux ratio. As can be seen, SN1997ey appears to have a cold SED, with a temperature of 20K or less, while it seems likely 

that SN1996cf's host is similar. 



Name 



RA 



DEC 



F450 mjy Fsso mJy L850 W m 2 sr 



SN1997cy 04 56 58.2 -02 37 37 0.58 20.8±3.5 7.8±1.1 
SN1996cf 10 48 50.9 00 03 30 0.57 25±15 11.0±1.6 
SN2000eh 04 15 02.4 04 23 18 0.49 5.2±6.6 4.6±1.5 



8.7xl0 22 
1.2xl0 23 
4.7xl0 22 



Table 2. Detailed Observational Properties of the strongly detected sources. For comnparison, Arp220, the nearest ULIRG, and Mrk231, 
the nearest type 1 AGN-ULIRG, have 850/im luminosities of 4.0 xl0 22 Wm~ 2 sr^ 1 and 1.0 X 10 23 tFm _2 sr _1 respectively. 



fluxes in the range 0.3 - 3 mjy at 850^im, and would thus be 
detectable by us either individually or statistically. And yet 
we find no evidence for such sources. Perhaps more signifi- 
cantly, sources as luminous as SN1996cf should be extremely 
rare - so rare in fact that we would not expect to see any 
at all given that we are only seeing 2 sources in the lumi- 



nosity range 10 22 ' 5 - 10 23 VKm _2 sr _1 . The assumption of a 
polynomial tail at high luminosities to the 850/im LF, rather 
than an exponential fall off, as suggested for the 60pm LF 
(Sanders & Mirabel, 1996) would ameliorate this problem 
to an extent, but would not provide a complete explaination 
for the numbers of high luminosity sources we are seeing, or 
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the absence of sources at somewhat lower luminosities that 
would still be detectable directly or statistically. 
Several possible explainations are possible: 

• The Luminosity Function has evolved from z=0 to 
z=0.5 

At face value, the most obvious explaination for the dis- 
crepencies in the luminosity distributions would be evolution 
in the luminosity function. This would not be one of the 
standard density or luminosity evolution since this would 
result in a many more moderate flux/luminosity detections 
than we see. It would also contradict the absence of evolu- 
tion in the opacity that we find for the sources undetected 
individually. Instead, some kind of evolution that only af- 
fects the high luminosity objects would be needed ie. some 
kind of luminosity-dependent luminosity or density evolu- 
tion. This could lead to a bimodal luminosity function. 

• The Local Luminosity Function is incorrect as a result 
of missing objects such as those we have detected 

While the SLUGS survey is currently the best determi- 
nation of the local 850pm luminosity function, it consists of 
less than 200 galaxies. It is thus entirely possible that classes 
of object are missing from it. If a class of cold dust temper- 
ature, high submm luminosity source existed, they might 
very well be absent from SLUGS. At the low dust tempera- 
tures of our galaxies, ~ 20K, such an object at z ~0.1 would 
have a 60/im flux ~ 70mJy, and so would not have appeared 
in the IRAS catalogs and would not have been targeted by 
SLUGS or any other submm followup. The corresponding 
100/im flux would be ~ 1 Jy, since this is reaching towards 
the peak of the SED. Such a cold source might be written 
off as likely cirrus contamination and thus be absent from 
the main IRAS galaxy catalogs. Its 850/im flux would be 
about lOOmJy, comparable to the 60/^m flux. At a redshift 
of 0.05 the fluxes would be 4 times larger and so such a 
source would still be only marginally detected by IRAS at 
60/im. If they exist, such objects would best be detected by 
large area, possibly all sky, shallow surveys with SCUBA2, 
or in the Planck all sky surveys. 

• The bright targets are not members of the same popu- 
lation as the rest of our targets 

Our application of the low z luminosity function is based 
on the assumption that all of our target SN host galax- 
ies come from the same parent population of established 
galaxies with large stellar populations. This might not be 
the case. If there was an additional population of SN host 
galaxies with large submm luminosities not arising from qui- 
escent dust heating, making up ~10% of SN hosts, then it 
would not be correct to extrapolate the submm properties 
of the bulk of the SN hosts from a few objects coming from 
a separate population. How might such a situation arise? 
The obvious class of objects with large submm luminosities 
are galaxies hosting a substantial burst of star formation, 
though the low dust temperature in our bright submm ob- 
jects and their quiescent morphology might argue against 
this. However, starburst galaxies are not the obvious hosts 
of type la supernovae. They are, though, the likely hosts of 
other classes of supernova. If some fraction of the SNla's 
in the surveys are actually misclassified core collapse super- 
novae in actively starforming, dusty hosts - they might per- 
haps be type lc SNe, which have masqueraded as type la's 
in the past (Kotak, private communication; see also Home- 



ier, 2004) - then this could provide an explanation. There 
would, though, be consequences for the SNla survey results 
if there is this level of contamination by non-la's. Core col- 
lapse supernovae are typically of lower luminosity than type 
las. Type lb and c SNe, viewed as a class, have absolute 
magnitudes about 1.24 magnitudes fainter than type la SNe 
(Richardson et al., 2002). If there is indeed a ~ 10% contam- 
ination rate of type las by type lcs in the surveys, then this 
could produce an overall dimming of the mean brightness 
by 0.24 magnitudes. This is comparable to the dimming at 
z~0.5 attributed to the cosmological constant or dark energy 
term (Perlmutter et al., 1999; Reiss et al., 1998), and thus 
might have a bearing on the reality or value of this cosmo- 
logical component. The redshift evolution in the dimming of 
SNlas attributed to dark energy at higher redshifts would 
not easily be mimicked by interloping type lc SNe, while the 
absence of starburst spectral features in these objects is an 
additional problem for this explanation. 

At this stage it is not possible to choose between these 
scenarios since the sample of such objects is still very small. 
Each possible explaination has its own significant implica- 
tions for astrophysics, so it is important that this strange 
population is followed up in more detail, and that steps are 
taken to expand the sample size. At this stage the only way 
to find more such objects seems to be to observe additional 
SNla host galaxies, in which they are represented at ~10%. 
This is not a particularly efficient method of discovery. Ex- 
tension of our SNla host studies to z=l might prove ben- 
eficial, and could test whether these objects are the result 
of evolution. Aside from that, the best hope for uncover- 
ing more of this class might be the cosmological submm 
surveys being proposed for SCUBA2. Furthermore, a very 
large area shallow SCUBA2 survey, possibly covering the 
whole accessible sky, would be one way to test the hypothe- 
sis that cold galaxies are a feature of the general galaxy pop- 
ulation at all redshifts by looking for their local equivalents. 
Conversely, detailed study of our existing submm bright tar- 
gets, in search of signs of recent star formation, might be the 
best way of testing our third explaination, and thus to re- 
move any lingering doubts about the classification of their 
supernovae. The absence of clear starburst indications or of 
morphological disturbance in existing data probably makes 
this final explanation the least favoured of the three we have 
proposed. 



6 CONCLUSIONS 

We have extended our search for dust in the host galaxies 
of type la supernovae at z~ 0.5. We combine our new ob- 
servations with those of Paper 1, and find that there is no 
evidence for a substantial increase in optical obscuration in 
the galaxy population between z=0 and z=0.5. We have, 
though, detected strong submm emission from one object, 
in addition to the one strong detection and one weaker de- 
tection from Paper 1. HST imaging of the two bright submm 
sources suggests they are quiescent, edge on disk galaxies, 
in contrast to the disturbed and interacting SMGs, found in 
blank field submm surveys, and local ULIRGs. We examine 
the dust properties of these sources, and conclude that they 
are ultra-luminous submm sources, with large dust masses, 
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but with cold dust SEDs, T~20K. Wo also examine the lumi- 
nosity distribution of the SNla host galaxies as a class, and 
compare it to the local SLUGS submm luminosity function. 
We conclude that the z—0.5 luminosity distribution seems 
bimodal, with an absence of sources within a factor of a few 
of our bright objects. Several scenarios could explain this ef- 
fect, including some forms of luminosity function evolution, 
the absence of sources such as these from local submm sur- 
veys, or problems with the classification of ~10% of SNe in 
the type la SN surveys. We propose various ways in which 
these hyptheses could be tested by future observations. 
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